Optical transitions and nature of Stokes shift in 
spherical CdS quantum dots 
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We study the structure of the energy spectra along with the character of the states participating in 
optical transitions in colloidal CdS quantum dots (QDs) using the ab initio accuracy charge patching 
method combined with the folded spectrum calculations of electronic structure of thousand-atom 
nanostructures. In particular, attention is paid to the nature of the large resonant Stokes shift 
observed in CdS quantum dots. We find that the top of the valence band state is bright, in contrast 
with the results of numerous k-p calculations, and determine the limits of applicability of the k-p 
approach. The calculated electron-hole exchange splitting suggests the spin-forbidden valence state 
may explain the nature of the "dark exciton" in CdS quantum dots. 

PACS numbers: 73.22.-f,71.15.Mb,79.60. Jv 
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Study of optical properties of semiconductor quantum 
dots is largely driven by potential applications afforded 
by their size-dependent bandgap and exciton spectra, 
such as solar cells^ lasers^ and fluorescent tags in by- 
otechnology applications. 3. In addition, quantum dots 
are excellent testing grounds for the applicability of vari- 
ous theoretical models. One of the most common features 
of quantum dots is the photoluminescence redshift rela- 
tive to absorption (also called Stokes shift). There are 
various causes for the Stokes shift, most common is due 
to non-resonant absorption in the existence of large size 
variation in the sample. However even when the non- 
resonant component has been eliminated, e.g, by the flo- 
rescence line narrowing techniques (FLN) where only the 
largest quantum dots are excited, there can still exist the 
Stokes shift. This resonant Stoke shift is usually caused 
by a dark exciton ground state, but the nature of this 
dark ground state can vary. A common dark ground state 
is caused by the exchange interaction between the elec- 
tron and the hole producing a spin triplet ground state 
which is spin-forbidden for optical transitions^ Another 
possibility consists of the electron and hole having dif- 
ferent spatial envelope function symmetries. Since it is 
relatively rare, this spatial symmetry induced dark exci- 
ton ground state has been actively pursued in quantum 
dots. It has been indicated in previous theoretical and 
experimental studies that the exciton ground state of a 
CdS quantum dot is one such candidatei^ In this work, 
we will use ab initio accuracy calculation to rc-investigate 
this problem. 

The use of the k-p theory to theoretically investigate 
the structure of energy levels and exciton states of spher- 
ical quantum dots has been very popularAiSiiiiiSiiiii It 
has been shown by various authors, including ourselves, 
tha t 7 i 8 i 14 based on the k-p calculations, the CdS QDs ex- 
citon ground state is an optically passive "dark exciton" 
for sufficiently small QDs because the hole ground state 
has an P-like envelope function. While the same have 
been found for CdSe and CdTe quantum dotsji^ the re- 
sults depend sensitively on the k-p parameters used. In 
contrast, for CdS quantum dots almost all reasonable k-p 



parameters, and even different k-p Hamiltonian models 
(including cubic and wurtzite crystal field splitting) pre- 
dict the same P-like hole ground state, thus indicating 
the spatial symmetry induced dark exciton. With in- 
creasing quantum dot size the hole ground state within 
the k-p framework eventually becomes a bright IS3/2 
state. This result has seemingly been confirmed by the 
large values of observed Stokes shift in CdS quantum 
dots; 7 ! 8 ! 15 ! 16 the calculated lP 3 / 2 -lS 3 /2 splitting values 
were shown to fit the observed values of resonant Stokes 
shift. 7 

While the k-p model had first successfully predicted 
the S-like state as the VBM in CdSe QDsi^, using a differ- 
ent set of k-p parameters generated from the pseudopo- 
tential bulk band structure it predicts the P-like state 
as the VBM, in contrast to the direct pseudopotential 
results 17 . This shows that the k-p and the pseudopoten- 
tial methods can predict different S/P ordering for the 
same set of k-p parameters. Unfortunately, the pseu- 
dopotential Hamiltonian for CdS quantum dot was pre- 
viously unavailable, therefore previous theoretical studies 
on CdS quantum dot have been limited to k-p-type cal- 
culations or tight-binding calculationsiSiiSiiSiiSii^ Re- 
cently, we have developed an ab initio accuracy charge 
patching method. This method can be used to calculate 
nanosystems of any given semiconductor material. A pre- 
vious quantum dot/quantum wire calculation using this 
method for thirteen different semiconductor materials^ 
have demonstrated an excellent agreement between the 
calculated and experimental optical band gaps. Here, we 
will apply this method to re-investigate the dark exciton 
problem in CdS quantum dots. 

The detailed recipe for our charge patching method 
calculations is given in Ref. Il9l here we only give a 
briefly outline. The CdS quantum dots are assumed to 
have wurtzite structure with lattice constants a = 4.12 
A and c = 6.73 A. The quantum dot effective radii 
R e f f is determined from estimating the number of atoms 
in a sphere of radius R e ff, assuming equal density of 
atoms in a quantum dot and the bulk. We are consider- 
ing quantum dots (CdS) 43 , (CdS)g 2 , (CdS)l83, (CdS)437, 



(CdS)874, (CdS)4586 (the subscript is the total number 
of Cd and S atoms), which have effective radii 6.32 A, 
8.15 A, 10.24 A, 13.70 A, 17.26 A, and 30 A respec- 
tively. The calculations were performed with the plane- 
wave pseudopotential method, using local density ap- 
proximation (LDA) and norm-conserving pseudopoten- 
tials. The plane wave energy cut-off of 35 Ry was used 
in all calculations. A well known LDA shortcoming of 
underestimating the bandgaps and the electron effective 
mass was corrected by modifying the nonlocal pseudopo- 
tentials (after Ref. Y19S) . such that the electron effective 
mass in the bulk is in a very good agreement with ex- 
periment, i.e. m e = 0.213mo calculated here, versus 
m e — 0.210mo measured^ and the bandgap is partially 
corrected (it is not possible to correct both simultane- 
ously), i.e. E g 



2.16 eV calculated versus E g = 2.58 eV 



measured** (uncorrected LDA yield; 
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0.127m 



and Eg — 1.315 eV). Spin-orbit coupling is included and 
adjusted to yield Ago=0.068 eV corresponding to the 
experimental valued In order to eliminate the quan- 
tum dot surface dangling bond states and keep the sys- 
tem charge neutral we passivate the surface with pseudo- 
hydrogen following Ref. 0. A surface atom of valency 
m is passivated with a pseudo-atom with Z — (8 — m) / 4, 
therefore Cd and S are passivated with H(Z = 1.5) and 
H(Z = 0.5), respectively. 

The charge patching method^ is used in order to ob- 
tain the self-consistent quality real space charge distribu- 
tion in the quantum dot, without having to perform di- 
rect LDA calculations, which for the sizes of the systems 
considered here are prohibitively expensive. Here, only 
small prototype systems are computed self-consistently 
for different atoms and their local environments to gen- 
erate the motif charge densities . These charge motifs 
are then used to assemble the total charge density for 
the entire quantum dot. Knowing the charge density the 
corresponding LDA potential is generated and the Hamil- 
tonian of a given quantum dot is constructed. The band 
edge eigenstates of this single particle Hamiltonian are 
then solved using the folded spectrum method^ 

Our calculated single particle eigenenergies for the va- 
lence band are shown in Fig^ Our calculation produces 
a IS3/2 state as the top of the valence band, and the 
energies in Fig^ are plotted relative to this state. For 
quantum dots in the region of Reff > 17A the IP3/2 is 
the second hole state. For smaller quantum dots, there 
is another S3/2 state above the IP3/2 state. Even for the 
largest QD calculated here, i.e. (CdS)4586 with effective 
radius of R e ff = 30A, we do not observe a S-to-P hole 
groundstate transition, and IS3/2 is still a groundstate. 
On the contrary, the 6x6 k-p theory in the spherical 
approximation, using the effective mass parameters de- 
rived from our ab initio bulk band structure (71 = 2.31, 
72 = 0.64, 73 = 0.89, A so = 0.068 eV, E p = 16.6 eV, 
Eg = 2.16 eV, m e = 0.21mo, and in the spherical approx- 
imation 72,3 = 0.79), yields the IP3/2 groundstate, rather 
than the lS 3 / 2 state, as shown in the inset to Fig. 2] For 
this set of parameters the P-to-S hole groundstate fran- 
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FIG. 1: Size dependence of the valence band energies without 
electron- hole Coulomb interaction in CdS quantum dots, plot- 
ted relative to the top of the valence band, IS3/2 state. The 
inset shows the hole energies of the IP3/2 and lS 3 / 2 states as a 
function of the quantum dot size, relative to the bulk valence 
band maximum, calculated with the spherical k-p method. 
The effective mass parameters used are, 71 = 2.31, 72 = 0.79, 
and A so = 0.068 eV. 



sition occurs at the QD effective radius of ^100A. This 
is similar to other k-p calculations even though k-p pa- 
rameters used might be somewhat different Ai2ii4 (Note, 
the caption of Fig. 5 in Ref. 7 is in error. The energy is 
given in the units of Eq = 7i/(2i? 2 ) not e = 7i/(2-R) 2 , 
and the x-axis is 0.529177 x i?(nm), not i?(nm).) Note, 
that the 8x8 k-p formula should give qualitatively sim- 
ilar results due to the large bandgap. Since the k-p pre- 
dicts IP3/2 hole groundstate for R e ff < 100A, our cal- 
culations indicate that it is possible that the size range 
of ~40A < R e ff < 100A, the IP3/2 is indeed the hole 
groundstate. However, the possible spatial symmetry in- 
duced dark exciton in that size range cannot be used to 
explain the experimentally observed Stokes shift, which is 
in the range of 15 - 70 meV, while the spatial dark exciton 
in that size range will have a Stokes shift of 0.6 - 6 meV, 
as predicted by the k-p model. The k-p method has been 
known to predict different ordering (and even omission) 
of the electronic levels in quantum dotsAZ, when compared 
to ab initio pseudopotential method. There has been a 
debate in the literature as to the proper comparison of 
the k-p and the ab initio methods 2 ^. While describing 
the bulk reasonably well, the two methods can give qual- 
itatively different results when applied to a nanocrystal 
(for a very detailed study see also Ref. |25J). Here we 
present another particular case of such disagreement. 

To compare our results with experimental optical mea- 
surements, in addition to the single particle energies, we 
need to calculate the exciton energies. Exciton energies 
of optical transitions in the strong confinement regime, 
where correlation effects are negligible, can be calculated 




FIG. 2: Size dependence of the energies of excitons con- 
structed from different valence band states including electron- 
hole Coulomb interaction. The inset shows the energy differ- 
ence between the IP3/2 and lS 3 / 2 valence state constructed 
excitons as a function of the CdS quantum dot size. 



FIG. 3: Electron-hole exchange splitting in comparison with 
experimentally measure resonant Stokes shift. Experimental 
data is deduced from Ref0 (Z. Yu et al.) and from Ref[ll (M. 
Chamarro et al.), measured with resonant photoluminescence 
(RPL) and photoluminescence excitation (PLE). 
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where, e c and e v are the single-particle conduction and 
valence states energies, respectively, and E^ v is the 
electron-hole Coulomb energy, obtained as2£ 
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where, x = (r, a) includes both spatial r and spin a =|, J, 
variables, e(ri — r 2 ) is a position-dependent dielectric 
function (described below), and V'c(x) and i/j v (x) are the 
wavefunctions for the conduction and valence states, re- 
spectively. The exciton energies can be further split by 
the electron- hole exchange interaction E* v in Eq.QJ. The 
exchange integral E* is calculated as: 
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In using the model dielectric function e(ri — r 2 ) we fol- 
low the procedure outlined in Ref. 1271 where the short 
range exchange interaction is essentially unscreened while 
the long range exchange interaction is screened signifi- 
cantly. Note that, to get the exchange splitting of the 
exciton energy, Eq.Q) needs to be diagonalized among a 
few Kramers doublet spin configurations, which are de- 
generate under the first three terms in Eq.JT]). 

The calculated exciton energies using Eq.Q without 
the E* v term are shown in Fig|21 The inset to FigEl 
shows the IS3/2 and IP3/2 energy difference which in- 
clude electron-hole Coulomb interaction as a function of 
the QD radius. The E\p 3 — £as 3/2 difference for the ex- 
citon has been enlarged compared to the single particle 



energies. This is because the lS 3 / 2 hole state has larger 
Coulomb interaction energy with the electron than the 
IP3/2 hole. Indeed, in the literature, there are reports 
of order changing in exciton states due to this Coulomb 
energy differenced In our case, it just makes the P-like 
hole state dark exciton less likely. 

We then add in the exchange interaction, and cal- 
culate the exciton exchange splitting for our calculated 
IS3/2 hole ground state exciton, therefore producing a 
spin-forbidden dark exciton. While there is no debate 
about the screening in the electron-hole Coulomb interac- 
tion, there exists a controversy regarding screening of the 
electron-hole exchange interaction in excitons. For ex- 
ample, when calculating optical properties in Si quantum 
dots 2 ^ and hydrogenated Si clusters 29 authors have calcu- 
lated exciton energy structure with unscreened exchange 
interaction, following the argument that in the Bcthc- 
Salpeter equation^ for the excitonic state the exchange 
term should be unscreened, otherwise leading to im- 
proper diagrams. Recently, however, it was shown^i that 
when the two-particle Green's function is constructed 
from a small set of particle-hole states, it may be appro- 
priate to screen the electron-hole exchange interaction. 
Therefore, in this work, along with the screened electron- 
hole exchange splitting we also calculate unscreened ex- 
change splitting for comparison. 

The results are summarized in Fig. [21 where we 
plot the singlet-triplet exciton energy splitting for the 
screened exchange interaction, unscreened exchange in- 
teraction, and experimental data taken from Refs. |7] and 
Il8l The overall trend of the electron-hole exchange split- 
ting agrees with the experimental data for the resonant 
Stokes shift. The experiment seems to favor the un- 
screened electron-hole exchange splitting results. This is 
in agreement with a previous tight-binding study«^ where 
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FIG. 4: Symmetric (a) and asymmetric (b) CdS quantum 
dots. The asymmetric dots exhibit a much (two orders of 
magnitude) smaller radiative lifetimes for P-type hole states, 
compared with symmetric dots. 

the exchange interaction was also unscreened. On the 
other hand, there is a systematic discrepancy between 
the experiment and calculated values of the screened 
electron-hole exchange splitting. Overall, our results in- 
dicate that singlet-triplet splitting could be responsible 
for the observed values of the resonant Stokes shift in 
the experiment. In order to confirm this theoretical pre- 
diction we propose an experiment in which the resonant 
Stokes shift is measured as a function of applied exter- 
nal magnetic field. The linear dependence of the Zeeman 
splitting on the external magnetic field could be detected 
from the Stokes shift measurement^ 

Last, to further study the possibility of the P-like hole 
induced spatial dark exciton, we have examined the ge- 
ometry dependence of the spatial dark exciton radiative 
life time. Note that the quantum dot (wurtzite) struc- 
tures considered so far are highly symmetric. Namely, 
there is a 120 degree rotation symmetry along the dot 
center z-axis (looking down in Fig. 0J. However, quan- 
tum dots obtained in experiment will always have some 
structural imperfections, away from the perfect symmet- 
ric shape. It is therefore interesting to study the influ- 
ence of such imperfection on the rate of optical transi- 
tions. Here we calculate radiative lifetimes r, following 
Ref. l32l Experimentally measured values for the radia- 
tive lifetimes are£ ~180 ns for the slow component of 



the luminescence at the temperature of 10 K which cor- 
responds to the radiative recombination of the optically 
passive state. In our calculations for the symmetric quan- 
tum dot [FigQJa)] we obtain radiative lifetimes of 1 ns for 
S'-like states, and very long (> 1000 ns) lifetime for P-like 
states, the latter is at least an order of magnitude larger 
than observed in experiment. We then calculated a sim- 
ilar sized quantum dot with its axial symmetry removed 
[Fig^Jb)]. As a result, the radiative lifetime of the dark 
IP3/2 induced exciton decreases sharply to ~20 ns, much 
smaller than the observed experimental lifetime. Since it 
is hard to expect perfect symmetric structures in the ex- 
periment, this result further supports the idea that the 
experimentally observed Stoke shift is not caused by the 
spatial dark exciton, rather it might be caused by the 
exchange interaction induced spin-forbidden dark exci- 
ton. Note that the spin-forbidden dark exciton will not 
be eliminated by the change of spatial shape of the quan- 
tum dot, and the CdS experimental lifetime of ^180 ns 
(QD of R = 11 A at T = 10 K) is similar to the dark 
exciton lifetime in CdSe (210 ns for DQ of R = 13 A at 
T = 12 K)2i where it is known that the dark exciton is 
caused by the exchange splitting. 

In summary, we have performed calculations of exci- 
ton states, optical properties, and electron-hole exchange 
splitting in CdS quantum dots. The results indicate that 
the previous k-p method wrongly assigned the top of va- 
lence band for small (R e ff < 30A) CdS quantum dots. 
Our hole ground state is found to be bright S-state, rather 
than dark P-state predicted by the k-p method. As a 
consequence, the dark exciton is not spatially forbidden 
in CdS quantum dot. Our calculation indicates that the 
exchange splitting might be responsible for the dark ex- 
citon and the observed Stoke shift. We recommend the 
magnetic field experiments to further resolve this issue. 
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